Intrinsic Magnetism in Nanosheets of Sn0 2 : A First-principles Study 
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We propose intrinsic magnetism in nanosheets of SnC^, based on first-principles calculations. 
The electronic structure and spin density reveal that p orbit als of the oxygen atoms, surrounding Sn 
vacancies, have a non itinerant nature which gives birth to localized magnetism. A giant decrease 
in defect formation energies of Sn vacancies in nanosheets is observed. We, therefore, believe that 
native defects can be stabilized without any chemical doping. Nanosheets of different thicknesses are 
also studied, and it is found that it is easier to create vacancies, which are magnetic, at the surface 
of the sheets. Sn02 nanosheets can, therefore, open new opportunities in the field of spintronics. 



The quest for room temperature (RT) ferromagnetism 
in diluted magnetic semiconductor (DMS), where tran- 
sition metal (TM) (magnetic) impurities are doped into 
semiconductor hosts, e.g., Mn in GaAs 1 , has motivated 
many theoretical and experimental studies. TM impu- 
rities were found however to produce intrinsic defects 
which limit its application in the field of spintronics at 
RT. Next ideal candidates for RT applications in spin- 
tronics are oxide-based DMS systems, which not only 
have high Tc, but also large magnetic moments 2 . Such 
oxide-based DMS also faced challenges, i.e., the deter- 
mination of the origin of their magnetic properties and 
the control of these properties in oxides. The difficulties 
arise from the presence of defects (cation/anion vacan- 
cies), which can be beneficial for magnetism but whose 
degree of influence is not completely understood. It is 
believed that some oxides show magnetism without any 
magnetic impuritieJ^' mainly due to cation vacancies. 
The role of vacancies (either cation or anion) in mag- 
netism can not then be ignorecP. 

For the last few years extensive work (both experimen- 
tal and theoretical) has been done to understand the ori- 
gin of magnetism in materials with defects. Some of these 
new discoveries also made possible to have magnetism in 
nonmagnetic insulators/semiconductors doped with light 
elements^UU. Very recently, some of us found that C can 
also induce magnetism in Sn02 9 , and later on this new 
theoretical idea was shown experimentally by Hong et 
afi2*. Therefore, it is now believed that magnetism can 
exist either in materials with vacancy or light elements- 
doped systems^H^Kl. Usually, cation vacancies, which 
are a source of magnetic moments, have large formation 
energies and it is not easy to produce them experimen- 
tally. Therefore, there remains an open and challenging 
question how intrinsic defects (vacancies) can be stabi- 
lized in nonmagnetic materials. Here, we find a new way 
of stabilizing defects in SnC>2. We show that intrinsic 
defects can be easily stabilized by reducing the dimen- 
sionality of Sn02 . This is a very unique way of reducing 
the formation energy of vacancies without any alien im- 
purities. Note that there are several experimental reports 
on the fabrication of nanostructured SnOsF^^. Beltran 



et a/P^, have discussed the thermodynamic stability of 
nanostructured Sn02 and have shown (both theoretically 
and experimentally) the possibility of nanosheets forma- 
tion. 

First-principles calculations based on density func- 
tional theory within the local spin density approxima- 
tion (LSDA/^ and generalized gradient approximation 
(GGA) 20 were performed using the SIESTA codepH 
which employs norm-conserving pseudopotentials and 
linear combinations of atomic orbitals. The convergence 
of the basis set (a double-zeta polarized, DZP) and other 
computational parameters (real space energy cutoff, 200 
Ry; fc-point sampling, 9x9x1 fc-points) were carefully 
checked. All atomic positions were fully relaxed until all 
atomic forces were smaller than 0.05 eV/A, which proved 
enough to converge the structure of the systemP Differ- 
ent single-layer nanosheets [2x2x1 (24 atoms), 3x3x1 
(54 atoms) and 4x4x1 (96 atoms)] were considered 
to take into account the effect of the vacancy concentra- 
tion on the electronic and magnetic properties (a 4 x 4 x 1 
nanosheet is shown in Fig.[TJa) and (b)). We also consid- 
ered nanosheets of different thicknesses [2 x 2 x n, where 
n = 2, 3, 4], as depicted in Fig. [I] to investigate the influ- 
ence of the position of Sn vacancy on the defect formation 
energies. The Sn and O vacancies were simulated within 
the supercell approach with periodic boundary conditions 
in the plane and open boundary conditions perpendicular 
to the plane; a vacuum region of 10 A along the direction 
normal to the surface of the nanosheet was found enough 
to achieve convergence on this parameter. 

First, we carried out nonmagnetic (NM) and mag- 
netic (M) calculations on pristine Sn02 nanosheets and 
we found no trace of magnetism. The total and atom- 
projected density of states (DOS and PDOS) are shown 
in Fig. [5|a). It is clear that the pristine sheet is a wide 
band gap nonmagnetic material where all bands below 
the Fermi level (2£p) are completely occupied, similar to 
bulk Sn02- The DOS near the Fermi level is dominated 
by the p orbitals of O in the valence band and the Sn 
orbitals (with some weight of the p O) in the conduc- 
tion band. Once confirmed that pure Sn02 nanosheets 
are nonmagnetic, we carried out simulations of native 
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FIG. 3. (Color online) Spin density of a single sheet SnC^ 
with a single Sn vacancy. Red (yellow) balls show the tin 
(oxygen) atoms. Arrows represent the direction of the spin 
magnetic moments. Blue (red) colors in the plane show neg- 
ative (positive) spin polarization of the atoms. 



FIG. 1. (Color online) Structure of nanosheets of Sn02- (a) 
shows a 4 x 4 x 1 nanosheet. (b) shows the side view of (a), 
(c), (d), and (e) show 2x2x2, 2x2x3, and 2x2x4 
nanosheets, respectively. Labels a, ...,<i indicate the location 
of the Sn vacancy. 



defects (i.e., Sn and O vacancies) in a single-layer Sn02 
nanosheets. We must note that in these nanosheets Sn 
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FIG. 2. (Color online)The total and projected density of 
states (PDOS) of pure (a) (left-side) and with a single Sn 
vacancy (b) (right-side) SnCb sheets. Solid and dashed lines 
show p and s orbitals of O/Sn atoms. Thick solid lines in 
the bottom panel show the total density of states. The Fermi 
energy (Ef) is set to zero. 

is not surrounded by six but by four oxygen atoms (see 
Fig. [TJb)). Then dangling bonds are expected in such 
structures. We created the Sn and O vacancies at the 
center of the sheets. Interestingly, we found that Sn va- 
cancies induce very large magnetic moments of 4.0/iB- 



However, the O vacancies do not induce magnetism, this 
behavior is similar to bulk SnO^. A detailed Mulliken 
analysis shows that the magnetic moment is mainly con- 
tributed by the oxygen atoms surrounding the Sn va- 
cancy. Each of the surrounding four O atoms carry a lo- 
cal magnetic moment of about 1.0/iB- This can be easily 
understood by taking into account that the Sn vacancy 
creates localized dangling bonds at the O sites which give 
rise to local magnetic moments following Hund's rules^. 
The magnetic coupling (ferromagnetic or antiferromag- 
netic) between the O atoms depends on the relative ori- 
entation and distance between the cation vacancie d 3 * 23 *. 
The localization of the O magnetic moment is further 
confirmed by LDOS and spin-density maps. 

Figure |2|b) shows the calculated total DOS and PDOS 
of a Sn02 sheet with a single Sn vacancy. We only show 
the PDOS of one of the oxygen atoms, since the remain- 
ing three oxygen atoms around the vacancy displayed 
the same behavior. The other oxygen atoms also con- 
tribute to the total DOS but their contributions are very 
small as compared with the oxygen atoms surrounding 
the Sn vacancy. The total DOS shows that the Sn va- 
cancy induces magnetism in Sn02 nanosheets but does 
not destroy the insulating behavior of Sn02- Therefore, 
Sn02 sheets with Sn vacancies can be assigned to the 
class of magnetic insulators. It is clear from this figure 
that the magnetization is mainly contributed by the oxy- 
gen atoms. The oxygen spin-up p band is completely 
occupied whereas the spin-down band has a smaller oc- 
cupation below Ey and some weight above it. The elec- 
tronic structure near the Fermi level is dominated by the 
oxygen p orbitals, but there is also a small contribution 
from the p and s orbitals of Sn below the Fermi level, 
where oxygen states hybridize with the Sn states. This 
partial/complete occupation of the oxygen p orbitals il- 
lustrates their localized nature. It must be noted that 
such localized p states are absent in bulk Sn02 3 , which 
suggest that this localized behavior of the p orbitals can 
be attributed to the reduced dimensionality of Sn02- 
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TABLE I. Calculated formation energies (in units of eV) in 
equilibrium (Eq.), Sn-rich, and O-rich conditions for a Sn va- 
cancy in Sn02 sheets as a function of Sn vacancy concen- 
tration (in %). The last column shows the total magnetic 
moment (MM) per unit cell (in units of //b)- 
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To shed more light on the nature of bonding and the 
source of magnetism in Sn02 nanosheets, we calculated 
the spin-density on a plane that contains the Sn va- 
cancy (see Fig. [3|. As can be seen, the spin-density is 
strongly localized at the oxygen atoms, which couple fer- 
romagnetically. The spin-density is, therefore, concen- 
trated around the Sn vacancy and does not extend to 
other O atoms, which signals the non-itinerant behav- 
ior of the oxygen p electrons. There are also small neg- 
ative induced magnetizations at the nearest Sn atomic 
sites. This shows that the negative polarization of the 
Sn atoms is essential for mediating magnetism between 
oxygens in the context of superexchange interactions^, 
which play a decisive role in governing localized mag- 
netism. We repeated the same calculations for different 
sheets (i.e., varying the O or Sn vacancy concentrations), 
and we found the same conclusion, i.e., the Sn vacancy 
induces magnetism, which is mainly contributed by the 
surrounding oxygen atoms. The magnetic moments are 
summarized in Table [H 

In order to determine the experimental feasibility of 
such configurations it is essential to determine the defect 
formation energy Ef. We calculated the formation ener- 
gies of the Sn and O vacancies in these nanosheets using 
the formula 

E { = £ D (Sn0 2 ) - £ P (Sn0 2 ) + ™Mx, 

where Ep (£b) is the total energy of pristine (defect, 
containing either Sn or O vacancy) sheets of Sn02 , /ix is 
the chemical potential of X(X = Sn or O), and n is num- 
ber of X atoms removed. In equilibrium conditions, the 
chemical potentials of O and Sn satisfy the relationship 
/iSn + 2/io = MSn0 2 ? where /isn0 2 > the chemical potential 
of bulk Sn02, is a constant value calculated as the total 
energy per Sn02 formula unit. Under Sn rich conditions, 
/i Sn =E(Sn metal ) and fi = ( /iSnQ2 2 ~ /fSn ), whereas in O 
rich conditions, ji = (^)E(0 2 ) and /i Sn = MSn02 - 2/io- 
The calculated Ef as a function of Sn vacancies, in 
equilibrium, Sn-rich and, O-rich conditions, is summa- 
rized in Table [IJ As can be seen, the formation of Sn 
vacancies is favored under O-rich conditions, and it de- 
creases as the concentration of Sn vacancies increases in 
the Sn02 nanosheets. It is very encouraging to see that 
Ef is significantly decreased in nanosheets as compared 
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FIG. 4. (Color online) Calculated defect (Vs n ) formation en- 
ergies Ef (eV) for Sn02 nanosheets with different thicknesses. 
The x-axis shows the location of Vs n in different sheets as 
marked in Fig [I] Square, circle, and triangle represent Ef 
of 2 x 2 x 4, 2 x 2 x 3, and 2 x 2 x 2 thick sheets of Sn0 2 , 
respectively. Filled (empty) symbols represent Ef in O-rich 
(Sn-rich) conditions. The horizontal dashed (solid) line shows 
Ef of bulk Sn0 2 in Sn-rich (O-rich) conditions. 



with Ef in bulk Sn0 2 , which is ~ 6.88 and 14.44 eV in 
O-rich and Sn-rich conditions, respectively. The smaller 
values of Ef in nanosheets of Sn02 (~ 75% smaller than 
the bulk values for the O-rich case) further guarantee the 
experimental realization of our proposed systems. We 
must also mention that previously synthesized samples of 
transition metal doped Sn02 have much larger Ef, e.g., 
7.0, 2.0 and 3.27 in O-rich conditions for Zn-, Cr-, and 
Co-doped Sn02, respectivel y 24 " 26 !. This further makes 
us confident about the experimental verification of our 
proposed systems. On the other hand, Ef of O vacan- 
cies in Sn02 nanosheets is calculated to be ~ —0.50 and 
~ 3.30 eV in Sn- and O-rich conditions, respectively. Al- 
though Ef is still smaller for O-vacancies in nanosheets as 
compared to the bulk, we believe these vacancies do not 
play a significant role in Sn02 in terms of d° magnetism, 
since they do not induce magnetic moments. However, 
we can not forget the role of O vacancies in promoting 
magnetism in TM-doped SnO^F 1 ^. 

Next, we considered nanosheets of different thicknesses 
[2x2x2,2x2x3, and 2x2x4] and created Sn vacancies 
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in different layers, as shown in Fig. [TJc)-(e). The calcu- 
lated formation energies are summarized in Fig. [4j The 
Sn vacancy created at the surface of the sheets [marked 
as a in Fig. [TJc)-(e)]) has the lowest formation energy, 
~ 0.70 eV, which is a giant decrease in Ef as compared 
with the bulk value. Such a small value further guaran- 
tees that Sn vacancies will be very active at the surface 
of the sheets and will easily be achieved experimentally. 
However, as the defect was moved to the second/third 
layer of the sheets [for example, sites b, c in Fig. [IJe)], 
E{ increased and approached the bulk value in both Sn- 
rich and O-rich conditions. This behavior can be un- 
derstood by taking into account bond saturation, which 
increases the stabilization of each atom inside the sheet. 
However, near the surface layers, the broken Sn-0 bonds 
reduce Ef significantly. In any case, whatever the lo- 
cation of the Sn vacancy in these sheets, Ef is always 
smaller than the bulk values. Once we have calculated 
Ef of the point defects, its concentration (c), which de- 
pends on its formation energy, in thermodynamic equi- 
librium is given by c = n s i te exp Where Ef is 

the defect formation energy, n s i te is the number of sites 
the defect can be incorporated on, &b is the Boltzmann 
constant, and T is the temperature. The above expres- 
sion shows that defects with high formation energies will 
occur in low concentrations. It is also very interesting 
that the magnetic moment (4.0/i#) remains robust when 
the thickness of the sheets changes or the Sn vacancy 
moves. We found, however, a significant increase in the 
local magnetic moments of the oxygen atoms at the sur- 
face of the sheets. The two surface O atoms have a local 
magnetic moment ~ 1.25/i#, which is almost indepen- 
dent of the sheet thickness. After relaxing the structure, 
the bond lengths of the Sn and O atoms surrounding the 
Sn vacancy were increased by ~ 0.07A. This relaxation 
was anisotropic, since the in-plane and out-of-plane bond 
lengths were slightly different. The structural relaxation 
had no significant effect on the total magnetic moments. 
However, a slight change in the local magnetic moments 
of the Sn and O atoms was observed. 

Our extensive calculations clearly demonstrate 
the possibility of magnetism without impurities in 
nanosheets of Sn02- This possibility is proved by 
the fact that the defect formation energies of the Sn 



vacancies are much smaller than those in bulk Sn02- 
This makes us believe that such magnetic nanosheet s are 
experimentally possible. Interestingly, Wang et a/ ! 27 * 28 * 
recently synthesized nanosheets of Sn02 and showed 
room temperature ferromagnetism in fresh and annealed 
systems. Such experimental report further confirms 
that room temperature ferromagnetism in nanosheets 
of Sn02 is enhanced by the small values of the defect 
formation energies. Finally, test calculations, using 
GGA and LDA+£7, were also carried out on pristine 
nanosheets and sheets with Sn vacancies. We found the 
same behavior, i.e., pristine nanosheets are nonmagnetic 
and Sn vacancies in nanosheets induce magnetism. Note 
that previous theoretical calculations have shown that 
the inclusion of the U term do es no t change the magnetic 
moments caused by vacancie d 29 * 30 *. Such term will have 
some effect on the electronic density of states^, but it 
will not change the main conclusion drawn here. Note 
that we focused on the origin of magnetism in Sn02 
nanosheets and o ur cal culations agree with the recent 
experimental dat d 27 * 28 -! The macroscopic magnetic 
coupling (FM or AFM) between Sn vacancies would 
required much larger simulations, which is beyond the 
scope of the present study. 

In summary, first-principles calculations based on den- 
sity functional theory within local density approximation 
have been carried out to investigate magnetism and for- 
mation energies of native defects in Sn02 nanosheets with 
different thicknesses. We showed that Sn vacancies in- 
duce large magnetic moments and O vacancies do not in- 
duce magnetism. The origin of localized magnetism was 
discussed using electronic densities of states and spin den- 
sity contours. The defect formation energies for different 
thicknesses were also calculated and it was found that Sn 
vacancies need much smaller energies in nanosheets than 
in bulk Sn02- These results indicate that it is easy to cre- 
ate Sn vacancies at the surfaces of Sn02 nanosheets. We 
conclude that intrinsics defects can be easily stabilized in 
nanosheets of Sn02 without any alien impurities. 

GR acknowledges the cluster facilities of NCP, Pak- 
istan. VMGS thanks the Spanish Ministerio de Ciencia 
e Innovacion for a Ramon y Cajal fellowship (RYC-2010- 
06053). JMM acknowledges the computational facilities 
provided by CESUP-UFRGS, Brazil. 



gulrahman@qau.edu.pk 

1 H. Ohno, A. Shen, F. Matsukura, A. Oiwa, A. Endo, S. 
Katsumoto, and Y. Iye, Appl. Phys. Lett. 69, 363 (1996). 

2 Ogale et a/., Phys. Rev. Lett. 91, 077205 (2003). 

3 G. Rahman, V. M. Garcia- Suarez, and S.C. Hong, Phys. 
Rev. B. 78, 184404 (2008). 

4 C. D. Pemmaraju and S. Sanvito, Phys. Rev. Lett. 94, 
217205 (2005). 

5 Q. Xu, H. Schmidt, S. Zhou, K. Potzger, M. Helm, H. 
Hochmuth, M. Lorenz, A. Setzer, P. Esquinazi, C. Mei- 
necke, and M. Grundmann,, J. Appl. Phys. 92, 082508 



(2008). 

6 M. Y. Ge, H. Wang, E. Z. Liu, J. F. Liu, J. Z. Jiang, Y. K. 
Li, Z. A. Xu, and H. Y. Li , Appl. Phys. Lett. 93, 062505 

(2008) . 

7 C. Moyses Araujo, Mukes Kapilashrami, Xu Jun, O. D. 
Jayakumar, Sandeep Nagar, Yan Wu, Cecilia Arhammar, 
Borje Johansson, Lyubov Belova, Rajeev Ahuja, Gillian 
A. Gehring, and K. V. Rao, Appl. Phys. Lett. 96, 232505 
(2010). 

8 Chang-wen Zhang and S-S Ya, J. Appl. Phys. 106, 063709 

(2009) . 



5 



9 G. Rahman and V. M. Garcia- Suarez, Appl. Phys. Lett. 
96, 052508 (2010). 

10 W.-Z. Xiao, L.-L. Wang, L. Xua, Q. Wan, and B. S. Zou, 
Solid State Commun. 149, 1304 (2009). 

11 C. Zhang, H. Kao, J. Dong, and P. Li, Phys. Status Solidi 
B 246, 1652 (2009). 

12 N. H. Hong, J.-H. Song, A. T. Raghavender, T. Asaeda, 
and M. Kurisu, Appl. Phys. Lett. 99, 052505 (2011). 

13 Z. R. Dai, Z. W. Pan, and Z. L. Wang, Solid State Com- 
mun. 118, 351 (2001). 

14 E. R. Leite, J. W. Gomes, M. M. Oliveira, E. J. H. Lee, 

E. Longo, J. A. Varela, C. A. Paskocimas, T. M. Boschi, 

F. Lanciotti, P. S. Pizani, and P. C. Soares, J. Nanosci. 
Nanotechnol. 2, 125 (2002). 

15 Z. R. Dai, J. L. Gole, J. D. Stout, and Z. L. Wang, J. Phys. 
Chem. B 106, 1274 (2002). 

16 M. Law, H. Kind, B. Messer, F. Kim, and P. Yang, Angew. 
Chem. 114, 2511 (2002). 

17 G. G. Khan, S. Ghosh, and K. Mandal, J. Solid State 
Chem. 186, 278 (2012). 

18 A. Beltran, J. Andres, E. Longo and E. R. Leite, Appl. 
Phys. Lett. 83, 635 (2003). 

19 D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566 
(1980). 

20 J. P. Perdew and K. Burke and M. Ernzerhof, Phys. Rev. 



Lett. 77, 3865 (1996). 

21 J. M. Soler, E. Artacho, J. Gale, A. Garcia, J. Junquera, 
P. Ordejon, and D. S. Portal, J. Phys.: Condens. Matter 
14, 2745 (2002). 

22 S. Blundell, Magnetism in condensed matter (Oxford Uni- 
versity Press, 2001) 

23 J. Osorio-Guillen, S. Lany, S.V. Barabash, and A. Zunger, 
Phys. Rev. Lett. 96, 107203 (2006). 

24 W. Wei, Y. Dai, M. Guo, K. Lai, and B. Huang, J. Appl. 
Phys. 108, 093901 (2010). 

25 W. Wei, Y. Dai, M. Guo, Z. Zhang, and B. Huang, J. Solid 
State Chemistry, 183, 3073 (2010). 

26 C. Zhang and S. Yan, J. Appl. Phys. 106, 063709 (2009). 

27 C. Wang, Mingyuan Ge, and J. Z. Jianga, Appl. Phys. 
Lett. 97, 042510 (2010). 

28 C. Wang, Q. Wu, H. L. Ge, T. Shang, and J. Z. Jiang, 
Nanotechnology 23, 075704 (2012). 

29 A. Droghetti, C. D. Pemmaraju, and S. Sanvito, Phys. 
Rev. B 78, 140404 R (2008). 

30 V. Fernandes, R. J. O. Mossanek, P. Schio, J. J. Klein, A. 
J. A. de Oliveira, W. A. Ortiz, N. Mattoso, J. Varalda, W. 
H. Schreiner, M. Abbate, and D. H. Mosca, Phys. Rev. B 
80, 035202 (2009). 

31 A. Droghetti, C. D. Pemmaraju, and S. Sanvito, Phys. 
Rev. B 81, 092403 (2010). 



